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According to the SchifF theorem the nuclear electric dipole moment (EDM) is screened in neutral 
atoms. In ions this screening is incomplete. We extend a derivation of the Schifl theorem to ions and 
molecules. The finite nuclear size effects are considered including Z^a^ corrections to the nuclear 
Schiff moment which are significant in all atoms and molecules of experimental interest. We show 
that in majority of ionized atoms the nuclear EDM contribution to the atomic EDM dominates 
while in molecules the contribution of the Schiff moment dominates. We also consider the screening 
of electron EDM in ions. 
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I. INTRODUCTION 

Permanent electric dipole moment (EDM) of elemen- 
tary particle or atom violates both P and T invariance. 
The Kobayashi-Maskawa mechanism leads to extremely 
small values of the EDMs of the particles. It is also too 
weak to explain the matter-antimatter asymmetry of the 
Universe. On the other hand, most of the popular ex- 
tensions predict much larger EDMs which are within ex- 
perimental reach. Therefore, measurements of EDM pro- 
vide an excellent method to search for physics beyond the 
Standard Model. The measurements of EDM in atomic 
and molecular experiments are presented in Refs. [ll-[l9|. 

The EDM of an atom is mostly due to either electron 
EDM and T,P-odd electron-nucleon interactions in para- 
magnetic systems (with non-zero electron angular mo- 
mentum J) or due to the T, P-odd nuclear forces in dia- 
magnetic systems (J = 0; nuclear-spin-dependent e-N 
interaction contributes here too) . The existence of T, P- 
odd nuclear forces leads to the T, P-odd nuclear moments 
in the expansion of the nuclear potential in powers of dis- 
tance R from the center of the nucleus. The lowest-order 
term in the expansion, the nuclear EDM, is unobservable 
in neutral atoms due to the total screening of the external 
electric field by atomic electrons [l^ . It might be possible 
however to observe the nuclear EDM in ions, where it is 
screened incompletely (see e.g. [2l| - |23j ). The first non- 
vanishing terms which survive the screening in neutral 
systems are the Schiff moment which was defined in Ref. 

(see also Refs. where the contribution of the 

proton EDM was considered) and the electric octupole 
moment (the latter vanishes in nuclei of experimental in- 
terest which have spin 1/2). More accurate treatment of 
the finite nuclear size in Ref. has shown that the 

atomic EDM is actually produced by the nuclear Local 
dipole moment which differs from the SchifF moment by 
a correction ~ Z^a^ where Z is the nuclear charge and a 
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is the fine structure constant. Since all experiments deal 
with heavy atoms this correction is significant. 

In the non-relativistic classical limit the screening for- 
mulas can be obtained in a very simple way. The second 
Newton law for the ion and its nucleus in the electric field 
reads 

(Mjv + N,m,)a, = (Z - N,)eEo (1) 
M^cln — ZcEn (2) 

mette = eEe, (3) 

where rrie and Mjv are the electron and nuclear masses; 
tti, aN and Oe are the ion, nucleus and electron average 
accelerations respectively, Eq is the external electric field. 
En is the average electric field at the nucleus, Ee is the 
average electric field at one of the ion electrons, e is the 
proton charge, N^. is the number of electrons in the ion. 
Since system of particles moves altogether, the averaged 
accelerations must be equal (a^ = a^y = Og), therefore 
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As we can see, the average electric field for electrons 
is suppressed by the ratio of masses rrie/MN that is 
very small for heavy atoms. It means that in the non- 
relativistic limit there is practically no effect related to 
the electron EDM in heavy atoms and ions, — dg ■ Eg w 0. 
The interaction of the nuclear EDM d with the external 
field, — d ■ En, is suppressed by the factor {Z — Ne)/Z. 

The same approach can be used to determine the elec- 
tric field at the nucleus in a diatomic molecule: 

(Ml -t- M2 + Neme)a, = {Zi + Z2 - Ne)eEa, 
M2a2 = Z2eE2N, 
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■Eq. (6) 
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Screening is stronger for diatomic molecules because of 
the factor M2/(A/i -1- A/2) that contains both nuclear 
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masses. This indicates that the nuclear motion can 
not be ignored. We also see that in neutral atoms and 
molecules the field at the nucleus is zero, therefore the in- 
teraction of the nuclear EDM d with the screened electric 
field vanishes, dEjv=0. Similarly, 



(Zi +Z2- Ne) 



Ml + M2 + Ne-me 



-En. 



(7) 



The different screening laws of EDM in neutral atoms, 
ions and molecules raise a number of new questions. For 
example, is the screening term in the nuclear Schiff mo- 
ment different in neutral atoms and ions? Can nuclear 
motion in molecules produce any additional effects which 
do not exist in a single atom? Are there any new effects of 
the electron density polarization in ions and molecules? 
Simple classical formulas presented above do not answer 
these questions. This motivates us to revisit the quan- 
tum Schiff theorem [2^ and extend it to the cases of ions 
and molecules. We also derive a formula which more 
accurately takes into account the finite nuclear size and 
calculate corrections to the nuclear Schiff moment. 

The present work is also motivated by new exper- 
iments. Effects of EDM in molecules are enhanced 
[25I [26I [28I [29I . This is why the molecular experiments 
are so popular now. Recently the EDM experiment has 
been started with molecular ions [l^. The EDM exper- 
iments with atomic ions in the storage rings have been 
considered too [191 . 



II. SCREENING OF EDM IN ATOMIC IONS 



Nuclear EDM and SchifF moment 



The charge distribution in a finite size nucleus can be 
written as p{r) = pq{v) -\- 5 p{t) , where J pod^r — 1, 6p{t) 
is due to the P, T-odd interactions. The P, T-odd term 
in charge density leads to the nonzero nuclear dipole mo- 
ment d = dl/I ~ Ze J d^rSpr, where Ze is the nucleus 
charge, e is the proton charge. Let us define N^, as the 
number of electrons. If A^e 7^ Z a system is an ion. In a 
neutral atom (N^ = Z) our derivation is expected to give 
the same results as the Schiff theorem [g^l including the 
effects of the finite nuclear size [13, [13 • 

The Hamiltonian of a single atom in an external elec- 
tric field Eq can be written in the following form: 



H^T + Vo + V + U + W, 



(8) 



where 



2me dK, 2 2Mn dqN ^ ' 

y = ^eR,Eo-ZeqA,Eo, 

dp{v) 



U = -Ze^ ^ j d- 



|R.j - qw - r| ' 



W = -dE, 
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Here and qAr are the radius-vectors of the electrons 
and nucleus correspondingly. The expression for U can 
be expanded in powers of r/Ri since the nuclei size is 
small compared to the atomic scales. Let us keep the 
first two nonvanishing terms: 

U ^-dey — -- 



In the above expansion the octupole term was omitted 
since it leads to the mixing of the states with high elec- 
tron angular momentum and its contribution to the total 
atomic EDM is smah [2i |. 

Following Schiff let us define the operator 

Q = i^- (9) 



Ze dq_ 



.N 



It is easy to check that there is a relation between [Q, Vq] 
and U 



U = 



Wo 



-47reS^V,(5(R, -qAr) (10) 

i 

S^^I^ZeJ d'Spr'r~^d J dVpoMr^j , (11) 

where the expression for the Schiff moment S has the 
same form as for a neutral atom ^241. Substituting ex- 
pression for U and W = Q,V\ into Eq. (j8]) we obtain 



H = Hn 



Q,H„ -47reS^V,;(5(R,;-qAr), (12) 



where Hq = T -|- Vq + ^ is the Hamiltonian of the system 
in the external electric field without P, T-odd terms. The 
calculation gives the following result for the commutator 
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Ze 
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where a^v is the nuclear acceleration operator. To obtain 
the average value of the acceleration operator we can use 
the Ehrenfest theorem: 



(F) ^ jZ - jVe)eEo 
Mn Mn 



(14) 



where F is the average force acting on the nucleus (see 
Eq. ([T])). Substituting the above expression to Eq. 
we obtain for the averaged commutator of Q and Hq the 
following equation 



1 



dEo. 



(15) 



Substituting this result into Eq. ([T^ we obtain the effec- 
tive Hamiltonian of the ion in the external electric field 



H = Hn 



f IdEo 



47reS V V,;(5(R, -qjv) 



(16) 

Note that the derivation above is done in the adiabatic 
approximation assuming that we can average over elec- 
tron motion when we calculate the nuclear motion, i.e. 
we assume irie <Si Mn- If the number of electrons Ne — Z 
the EDM term in the above expression vanishes, as the 
SchifF theorem predicts. In the ion case the nuclear EDM 
interacts with the average field Ej^ = Eq{1 — Ne/Z) that 
acts on the ion's nucleus. 
The last term in Eq. (fTH]) . 



= -47reS ^ Vi<5(Ri - q^), 



(17) 



induces the ion EDM directed along the nuclear spin 
(which is the direction of the nuclear Schiff moment S), 
similar to the EDM of neutral atoms. This expression 
is not applicable for heavy atoms where the Dirac equa- 
tion gives infinite results for the electron wave functions 
at the point-like nucleus. Accurate account of the fi- 
nite nuclear size gives the following form for the cor- 
rected Schiff moment electrostatic potential (defined by 
= -e(/35(R)): 



3S' R 



B 



(18) 



where B = j po{R)R'^dR is the normalization constant. 
In the limit of the point-like nucleus the expression (jl8p 
agrees with Eq. pT)) . The corrected Schiff moment S' is 
given by the equation (see Appendix) 



S' 



Ze 
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(19) 



(rO--(r)(r^)--(n)(g,,r2) 



where qij is the quadrupole moment tensor. Here we 
omitted higher order terms which are proportional to a 
small factor Z'^a'^/9. Outside the nuclear radius Rn the 
nuclear density po (R) — and the potential ([TS]) vanishes 
in agreement with the Schiff theorem. Near the origin 
Po (R) — const and the potential (jlSp is a linear function 
of R. Therefore, the gradient of the Schiff moment poten- 
tial (IT51) gives a constant electric field inside the nucleus 
which is directed along the nuclear spin. This electric 
field polarizes the electron distribution and produces the 
atomic EDM. The calculations of the atomic EDM have 
been performed, for example, in Refs. (23. f30l - l32j . 

Below we make rough estimates to compare the nuclear 
EDM and the Schiff moment contributions to the atomic 
EDM. In the case of a spherical nucleus the nuclear EDM 
d, the nuclear Schiff moment and the atomic EDM Da 
induced by the Schiff moment have been estimated in 
Ref. M: 



10 ^^rje ■ cm , 



Da ^ {Z/lOOf ■ IQ-^'^rje ■ cm , 



(20) 
(21) 

where r] is the strength constant of the nuclear P, T-odd 
interaction (in units of the weak Fermi constant G) . As- 
suming the single ionization we get for the nuclear EDM 
screening factor 1 — Ng/Z = 1/Z. As a result, for the 
ionic EDM induced by the nuclear EDM we get the es- 
timate 1/Z ■ 10~^^77|e|cm. Thus, for the spherical nu- 
clei the nuclear EDM contribution exceeds the nuclear 
Schiff moment contribution by at least one order of mag- 
nitude. However, in heavy ions containing nuclei with 
the octupole deformation (e.g. ^^^Ra"*" and ^^'^Rn+) the 
Schiff moment contribution is enhanced by three orders 
of magnitude jssl . [33 | and is comparable to the nuclear 
EDM contribution (which is also slightly enhanced in 
these ions). 



B. Electron EDM 

For neutral atoms the electron EDM problem was in- 
vestigated in [sl] and further developed in [s!]. The 
Hamiltonian of the nucleus and relativistic electrons in 
the external electric field £^0 can be presented as 



Ho = -h'^AN/2MN - ZeqAfEo- 



(22) 
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-ihcaiVi + A 
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(23) 
(24) 



S = 
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where Et is the total electric field acting on the electron 
which includes the external field Eg, the nuclear field 
and the field of other electrons, a and /3 are the Dirac 
matrices. It is convenient to present as the sum of 
two terms 



Hui — Hid + , 
Hid = -de £ S,Et 



H 



2d 



AT, 



(A - l)S,Et 



(25) 
(26) 

(27) 



As it was pointed in [35|| the first term Hid gives no 
contribution to atomic EDM in a neutral atom. In an 
ion the Hid contribution is suppressed by a small factor 
me/Miq. It can be shown using the commutator relation 



Hid — Q, Hq 



We 



e ^ oKi 



(28) 
(29) 



Note that the matrix elements of the operators in the 
the Hid come from the atomic size area where valence 
electrons (which contribute to the atomic angular mo- 
mentum and EDM) are non-relativistic. To estimate the 



average value of the commutator Q, Hq 
theorem can be employed 



the Erehnfest 



[[Q,Ho 



(30) 



Substituting expression ([S]) for Eg into above equation 
we obtain for the average value of Hid 

{Hid) ^ -de^{Z - Ne){Y,^^^^) (31) 

i 

We see that the averaged value {Hid) is suppressed by 
the small mass ratio me/Mjq. It means, that in the limit 
of heavy nucleus Hid gives no coiitribution to EDM. 

The second perturbation term H2d vanishes in the non- 
relativistic limit since the matrix {Pi — l) acts on the lower 
components of the Dirac 4-spinors only. The operator 
H2d induces atomic EDM given by the same expression 
as for neutral atoms, except for the sum in the matrix 
elements is taken over electron number < Z: 

D2 =4(0|^(/3.-l)i:.|0) + 

E(ft-l)S.E,|n)(n|ER.|0) (32) 



2ede ^ ■ 



En — En 



In heavy atoms the major contribution to D2 comes from 
the second term {D2 ^ 3RreiZ^a^de where Rrei ~ 3 is 



the relativistic factor [35|,l36|) . This term corresponds to 
the atomic EDM due to the perturbation of the electron 
density by the operator H2d- Note that a similar equation 
with the perturbation Hid gives zero result due to exact 
cancellation between the first and second terms. Indeed, 
the zero and the first order corrections to the atomic 
EDM induced by Hid give EDM 

Di =de(0|^S,|0) + 

(0| [q,Ho\ |n)(n|ER»|0) 



^E- 

n 

=E 



En — En 



(0|ER»("l Q,Ho 



(33) 



|0) 



Eq — En 

The above expression can be simplified in the following 
way. For the matrix elements of the commutators the 
following relations are valid 



g,Ho \0) = -{Eo-En){n\Q\0) 



(0| 



Q: Hq 



\n) = {Eo - En){0\Q\n) 



(34) 
(35) 



Substituting these expressions into Eq. ([5^ and using 
the completeness condition ^ = 1 we obtain 



Di =e^^ [(0|g|n)(n|R,|0)-(0|R,|n)(n|Q|0)' 



+ 



|^S,;|0) =4(0|S,;|0)+^e( 



|o) 



(36) 



Using definition of the operator Q it is easy to show that 

QjRi = —de/e'Si. Hence, the second term in the above 

equation cancels the first term, so the dipole moment Di 
induced by Hid equals to zero. In this derivation we as- 
sume that the electron states are stationary. This is valid 
if we neglect the ion acceleration. Therefore, the result 
is consistent with Eq. (PT|) . 

We see that EDM of an ion induced by the electron 
EDM is given by the same equation as for neutral 
atoms (up to corrections ^ me/Mjv). A similar conclu- 
sion is also valid for molecular ions. 



III. NUCLEAR EDM AND SCHIFF MOMENT 
IN MOLECULAR IONS 



Let us consider a molecular ion with electrons and 
two nuclei with charges Zie and Z2e. We assume that 
the second nucleus has EDM d and Schiff moment S. 
The molecular Hamiltonian is equal to the sum of the 



5 



following terms 

We 



2me 2 2Mi ^qi 2 21/3 9q2 ^ ' 
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V = ^eR,;Eo - ZieqiEo - Z2eq2Eo, 



JVe 
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- q2 - rl' 



+ ZiZ2e' j d^r 
W = -dEo, 



2 / 



|qi - q2 - r| 



where qi and q2 are the coordinates of first and second 
nuclei respectively. Using the operator 



Z2e dq2 



(37) 



we can present the molecular Hamiltonian in the form 
similar to Eq. (fT2|) : 



(38) 



We 



47res <j - q2) - ^1 - q2) 



To calculate the average value of the commutator Q and 
Hq we can use the same algorithm as for a single atom. 



d i 

h 



= -^M2k2 (39) 



Since the molecule moves as a single body the average 
accelerations of all its particles is equal to the molecular 
acceleration, i.e. 



(a2) 



(F) 



(Zi + Z2- iVe)eEo 



Ml + M2 + NeTUe Ml + M2 

M2 Z1 + Z2- N, 



(40) 



Q, Hq 



dEo. (41) 



Ml + M2 Z2 
Finally, the effective Hamiltonian of the molecular ion is 
M2 Z1 + Z2- N, 



H=Hn- 



Mi + M2 

We 



-dEo 



(42) 
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AneS \ V,<5(R,; - qz) - Zig^^Sim - q2 



Thus, in a molecular ion the EDM term experiences 
the extra suppression. As for the Schiff moment term, 
it is still described by the same operator as for a sin- 
gle atom, except for the extra term proportional to 
d{S{qi —(I2))/ 9qi describing the interaction of the charge 
of the first nucleus and the Schiff moment of the second 
nucleus. The matrix elements of such interaction are ex- 
tremely small due to the Coulomb barrier. 



IV. ENHANCEMENT OF THE SCHIFF 
MOMENT CONTRIBUTION TO P,r-ODD 
EFFECTS IN POLAR MOLECULES 

Now we can compare the contributions of the nuclear 
EDM and Schiff moment to P, T-odd effects in polar 
molecular ions. Important difference between molecules 
and single atoms is that the nuclear motion significantly 
affects induced P, T-odd effects. The Schiff moment con- 
tribution in polar molecules is enhanced because of the 
strong internal electric field |25|]. Another interpretation 
of the enhancement is due to the small distance between 
the opposite parity rotational levels [l^, [3 ■ 

The nuclear P, T-odd effects are studied in the 
molecules with zero electron angular momentum. After 
averaging Hamiltonian Eq. (1421) over electron wave func- 
tion we obtain the effective Hamiltonian for the nuclear 
motion: 



2/i 



{q-qef+BJ{J+l) + Hu, (43) 



where q = qi — q2 , is the equilibrium distance between 
the nuclei in averaged potential, J is the rotational an- 
gular momentum of the molecule, Ue describes the inter- 
action of the partially screened nuclear EDM, the Schiff 
moment term H^j can be presented as [13, HI] 



H,„ =6XS 



In 



(44) 



where S = SI/ 1, n is the unit vector along the molecular 
axis, X is the constant that appears after averaging the 
perturbation over the electron wave function. In the first 
order of the perturbation theory the Schiff term leads to 
the rotation state mixing 



I ^ 



Ej - Ej 



(45) 



where ■0'-'^) — \Jm) is the unperturbed rota- 
tional wave function. Since the energy difference 
Ej-Ej, = B{J{J+\)-J'{J' + l)} can be very smah for 
rotation levels, the state mixing can be significant. This 
mixing induces EDM in the rotational state 



QXSDmI 



J(J + 1) - 3m2 



IB 



J(J + 1)(2J 



l)(2J + 3) 



(46) 



(47) 
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Here Dm = DmJi is the internal EDM of the polar 
molecule. This formula is valid for J ^ 0. For J = 
the induced EDM is 



2XSDmI: 

Tb 



(48) 



There is also the screened nuclear EDM contribution 
to P, T-odd molecular EDM ( see Eq. ^ Combin- 
ing this contribution with the Schiff moment contribution 
Z)f we obtained the P, T-odd part of the interaction of a 
molecular ion with the external electric field Eq: 



V = - 



AhTW, 'i-^-S — (49) 



This equation tells us that there is actually no enhance- 
ment of the electric field in the polar molecule since 
the electric field at the nucleus is suppressed I/Z2 times 
rather than enhanced. However, there is huge enhance- 
ment of the Schiff moment contribution since the expres- 
sion for the coefficient Km contains in the denominator 
the rotational constant B which may be five orders of 
magnitude smaller than the interval between atomic lev- 
els of opposite parity. 

Note that we can derive Eq. p9|) treating Eq as a 
perturbation. Therefore, the energy shift produced by 
the Schiff moment in Eq. is actually proportional 

to the average polarization of the polar molecule in the 
electric field Eq. In the small electric field it is linear in 
Eq, however, in the high field it tends to the constant. 
This determines the saturation effect in the energy shift 
produced by the Schiff moment if we go beyond the weak 
electric field Eq approximation (see Eq. (|44| where the 
average polarization < 1) . 

Using Eq. we can compare molecular EDM in- 

duced by the screened nuclear EDM and the Schiff mo- 
ment. Consider, for example, molecule PbF'^ since it has 
the same number of electrons as a well studied molecule 
TIF where the effect of the nuclear Schiff moment has 
been measured. The screened EDM term for PbF+ is 
Dn ^ 10-23,ye-cm ( EDM of F and EDM of odd iso- 
tope of Pb give comparable contributions since values of 
M/Z are approximately the same). To obtain the Schiff 
moment induced EDM in the ground state we need to 
estimate the constant K^, given by Eq. (^5)) . Since 
the molecular parameters are unknown for the ion we 
assume them to be of the order of their values for the 
neutral molecule TIF: X « 8000 a.u. [H, [13 . The 
values of the rotational constant B = 1.025 • 10~^a.u. 
and dipole moment Dm = 1.65 a.u. for TIF are taken 
from 38]. Finally, substituting all the parameters into 
Eq. (|48)l we obtain Km = 5 ■ 10^° a.u. Assuming the 
Schiff moment value for an odd isotope of Pb equal to 
S — 10~*77e • fm^ [111 we obtain the value for the Schiff 
moment contribution Ds ^ lO^^'^rje ■ cm which is three 
orders of magnitude larger than the nuclear EDM contri- 
bution Dm ^ 10"^'^ rye ■ cm. As it was mentioned above, 
in the nuclei with the octupole deformation like Ra^^^ 
the Schiff moment is enhanced. Therefore, in molecular 



ions like RaF"*" the Schiff moment induced EDM will be 
5 orders of magnitude larger than the partially screened 
nuclear EDM. 



V. CONCLUSIONS 

Accurate treatment of the electron EDM effects shows 
that the T,P-odd EDM of atomic and molecular ions at 
high Z are dominated by the Z^ enhanced relativistic 
correction effect, similar to neutral systems. The direct 
contribution of electron EDM is suppressed by the screen- 
ing factor {rrie/M) where M is the ion mass. 

The situation is different for the nuclear EDM. In 
atoms the nuclear EDM is screened by the factor Zi/Z 
where Zi is the ion charge. However, the nuclear EDM 
still dominates over the Schiff moment induced atomic 
EDM (with exception of heavy ions which contain nu- 
clei with the octupole deformation like ^^^Ra and ^^"^Rn 
where the Schiff moment is strongly enhanced). 

In molecular ions the nuclear EDM screening is slightly 
stronger than in atomic ions, the screening factor is 
{Mn /AI){Zi/Z). At the same the Schiff moment con- 
tribution is enhanced ~ Misf/me ^ 10^ times due to the 
mixing of the close rotational states of opposite parity. 
There is the additional Schiff moment enhancement in 
such molecular ions like RaF+. As a result, the Schiff 
moment contribution is 10^ — 10^ times larger than the 
screened nuclear EDM contribution. 

This combination of the large enhancement factors 
makes molecular ion experiments an attractive alterna- 
tive to the atomic EDM experiments. 



VI. APPENDIX 

According to Eq. ^TE\i in the limit of the point-like nu- 
cleus the Schiff moment potential and its matrix element 
are given by 

(^s(R) =47rS- V(5(R) (50) 
(s| - e^s\p) = 4^eS • (Vi^t^p)fl^=o (51) 

For the solutions of the Dirac equation (V-!/'JV'p)fl:->-o is 
infinite for a point-like nucleus. Therefore, for relativis- 
tic electrons it is necessary to account for the finite size 
of the nucleus and introduce a finite-size Schiff moment 
potential. An appropriate potential has been shown [27| 
to increase linearly inside the nucleus and vanish at the 
nuclear surface: 



3S' ■ R 

B 



7l(i?), 



(52) 



where B = J n{R)R'^dR « R%/5, Rn is the nuclear 
radius and n(R) is a smooth function which is 1 for 
R < Rn — S and for i? > R^ + S; n[R) can be taken as 
proportional to the nuclear density po (note that we can 
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choose any normalization of n(r) since the normalization 
constant cancels out in the ratio n/B, see Eq. ([5^ ). 

Below we will accurately derive expression for the cor- 
rected Schiff moment S' that corresponds to the potential 

The P, T-odd part of the nuclear electrostatic potential 
with electron screening taken into account can be written 
in the following form (see e.g. [sJl for the derivation): 



ep(r) 
R-r 



d^r (53) 



As it was shown in [27| the expansion of the Coulomb 
potential in (j53p in terms of the Legendre polynomials 
gives the following dipole term in the potential: 



i?3 



r 



r 

l3 



(54) 

We see that (^(^^(R) = if i? > Rn (nuclear radius) 
since /3(R) = in that region. Therefore, correspond- 
ing matrix elements will depend on the electron wave 
functions behavior inside the nucleus. All the electron 
orbitals for I > 1 are extremely small inside the nucleus. 
Therefore, we can limit our consideration to the matrix 
elements between s and p Dirac orbitals. We will use the 
following notations for the electron wavef unctions: 



^(R) = 



Lm 

-i{<T ■ n)g{R)n 



(55) 



where fJj/m is a spherical spinor, n = R/i?, f{R) and 
g{R) are the radial functions. Using (cr • n)^ = 1 we can 
write the electron transition density as 

PspiR) = i^li^p = nliipUspiR) (56) 

oo 

Usp{R) ^fs{R)fp{R) + 9s{R)9p{R) = J2 ^f^^" (57) 

The expansion coefficients bk can be calculated analyti- 
cally |27i]; the summation is carried over odd powers of 
k. Using Eqs. (I54I56|) we can find the matrix elements 
of the electron-nucleus interaction, 



{s\ - e(p(i)(R)|p) = ~Ze'{s\n\p) ■ |^°° [((r) - r) 
, {n)qij 



UspdR 



Ze'^{s\n\p) 



UspR dR 




k + 4: 



(58) 



where {s\n\p) = J Q,lni}pd(l) sin 0d9, (r") = / p{r)r"d^r. 
Note, that all vector values (rr") are due to P, T-odd 
correction Sp to the nuclear charge density po, while (r") 
are the usual P, T-even moments of the charge density 
starting from the mean-square radius (r^) = for k — 1. 

We now set the matrix elements ([58| of the true nuclear 
T, P-odd potential to be equal to the matrix elements of 
the equivalent potential (|5^ which are given by 



e0(R)|p) = 15e(s|n|p) • / UspR^n{R)dR 



15e(s|n|p) • 



R 



fe-i 

N 



fe=l 



fc-l-4' 



(59) 



where we 
/ n{R)R''dR P 
(15^1) we obtain 



Ze 



have 
/{k 



made approximation 
4). Equating dSH]) and 



15 ^^7? 



1 



k=l 



6i fc + 1 



(rr'^+i) - (r)(r''+^) 



k+l\ 



1 



. (60) 



Thus we have a possibility of separating the nuclear 
and electronic parts of the calculation of atomic EDMs. 
The nuclear calculation involves only the determination 
of S' and the atomic calculation involves only the effects 
produced by the equivalent potential ([5^ . 

Note that S' in eq. (pO]) is different from the Local 
dipole moment L defined in Ref. [131 : L does not contain 
the sum in the denominator. The reason for the difference 
is that here we reduce the problem to the nuclear size 
effective potential (|52l) while in Ref. 27[ the problem was 
reduced to the contact effective potential (|50p located in 
the center of the nucleus. 



In the non-relativistic case {Za 
hi ^ 0, and 



0) we have just 



hm S' = ^ 

Za^O 10 



{^r^)-\{^){r^)-\{n){m) 



(61) 



This is the usual expression for the Schiff moment S. In 
practice it may be sufficient to use only the first and 
third terms in the expansion of Ugp- In this case we 
need only the ratio 63/61. This ratio is different for the 
matrix elements s - pi/2 (63/61 = — (3/5)Z^Q;^/i?^) and 
s - P3/2 (63/61 = -{9/2Q)Z'^a^/R%). However, with the 
10% accuracy we can use the average of these two values 
63/61 « -Q.bZ^a^/R%. 



[1] D. DeMille et al, Phys. Rev. A. 61, 052507 (2000). 



[2] L. R. Hunter et al., Phys. Rev. 132, 2194 (2002). 



8 



[3] D. Kawall, F. Bay, S. Bickman, Y. Jiang, and D. DeMille, 

Phys. Rev. Lett.' 92, 133007 (2004). 
[4] M. G. Kozlov and D. DeMille, Phys. Rev. Lett. 89, 

133001 (2002). 

[5] A. N. Petrov, A. V. Titov, T. A. Isaev, N. S. Mosyagin, 

and D. DeMille, Phys. Rev. 132, 2194 (1963). 
[6] J. J. Hudson, B. E. Sauer, M. R. Tarbutt, and 

E. A. Hinds, Phys. Rev. Lett. 89, 023003 (2002). 
[7] B. E. Sauer, H. T. Ashworth, J. J. Hudson, M. R. Tarbutt 

and E. A. Hinds, Atomic Physics 20, 44 (2006). 
[8] N. E. Shafer-Ray, Phys. Rev. A 73, 034102 (2006). 
[9] M. G. Kozlov,V. L Fomichev, Y. Y. Dmitriev, L. N. Lab- 
zovsky, and A. V. Titov, J. Phys. B. 20, 4939 (1987). 
[10] Y. Y. Dmitriev at al, Phys. Lett. A 167, 280 (1992). 
[11] A. C. Vutha at al, J. Phys. B 43, 074007 (2010). 
[12] E. R. Meyer and J. L. Bohn, Phys. Rev. A 78, 010502(R) 
(2008). 

[13] J. Lee, E. R. Meyer, R. Paudel, J. L. Bohn, and 
A. E. Leanhardt, J. Mod. Opt. 56, 2005 (2009). 

[14] B. C. Regan, E. D. Commins, C. J. Schmidt, and D. De- 
Mille, Phys. Rev. Lett. 88, 071805 (2002). 

[15] J. M. Amini, C. T. Munger, and H. Gould, Phys. Rev. 
A 75, 063416 (2007). 

[16] F. Fang and D. S. Weiss, Opt. Lett. 34, 169 (2009). 

[17] S. A. Murthy, D. Krause Jr., Z. L. Li, and L. R. Hunter, 
Phys. Rev. Lett. 63, 965 (1989). 

[18] A. E. Leandhardt, lar^v :1008.2997 (2011). 

[19] Workshop on search for electric dipole moments (EDMs) 
at storage rings, 

http: / / www2.fz-juelich.de / ikp /edm / en /program. shtml 
[20] L L Schiff, Phys. Rev. 132, 2194 (1963). 
[21] V. A. Dzuba, V. V. Flambaum, P. G. Silvestrov and 

O. P. Sushkov, Phys. Lett. A 118(4), 177 (1986). 
[22] C. Itoi, S. Oshima, The Phys. Soc. of Japan 79(10), 



103201 (2010). 

[23] V. A. Dzuba, V. V. Flambaum, 'arXiv:1009.4960l(2010). 
[24] O. P. Sushkov, V. V. Flambaum, I. B. Khriplovich, Zh. 

Ex. Teor. Fiz. 87, 1521 (1984) [JETP 60, 873 (1984)]; 

V. V. Flambaum, I. B. Khriplovich, and O. P. Sushkov, 

Nucl. Phys. A 449, 750 (1986). 
[25] P. G. H. Sandars, Phys. Rev. Lett. 19(24), 1396 (1967). 
[26] E. A. Hinds, P. G. H. Sandars, Phys. Rev. A 21(2), 471 

(1980). 

[27] V. V. Flambaum, J. S. M. Ginges, Phys. Rev. A 65, 
032113 (2002). 

[28] D. Cho, K. Sangster, E. A. Hinds, Phys. Rev. A 44(5), 
2783 (1991). 

[29] O.P. Sushkov, V.V.Flambaum. Zh.Exp.Teor.Fiz 75, 1208 

(1978) [Sov.Phys. JETP 48, 608 (1978)]. 
[30] V. A. Dzuba, V. V. Flambaum, J. S. M. Ginges, Phys. 

Rev. A 61, 062509 (2000). 
[31] S. G. Porsev, J. S. M. Ginges, V. V. Flambaum, Phys. 

Rev. A 83, 042507 (2011). 
[32] V. A. Dzuba, V. V. Flambaum, J. S. M. Ginges, 

M. G. Kozlov, Phys. Rev. A 66, 012111 (2002). 
[33] N. Auerbach, V. V. Flambaum and V. Spevak, Phys. 

Rev. Lett. 76, 4316 (1996). 
[34] V. Spevak, N. Auerbach, V. V. Flambaum, Phys. Rev. 

C 56, 1357 (1997). 
[35] P. G. H. Sandars, Phys. Lett. 14, 194 (1965). 
[36] V. V. Flambaum, Yad. Fiz. 24, 383 (1976) [Sov. J. Nucl. 

Phys. 24, 199 (1976)]. 
[37] A. N. Petrov et al, Phys. Rev. Lett. 88, 073001 (2002). 
[38] K. P. Huber, G. Herzberg, Molecular Spectra and Molec- 
ular Structure IV. Constants of Diatomic Molecules (Van 

Nostrand: Reinhold, New York, 1979). 



